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We present a pressure sensor based on a Michelson interferometer, for use in photoacoustic tomography. 
Quadrature phase detection is employed allowing measurement at any point on the mirror surface without 
having to retune the interferometer, as is typically required by Fabry- Perot type detectors. This opens the door 
to rapid full surface detection, which is necessary for clinical applications. Theory relating acoustic pressure 
to detected acoustic particle displacements is used to calculate the detector sensitivity, which is validated 
with measurement. Proof-of-concept tomographic images of blood vessel phantoms have been taken with 
sub-millimeter resolution at depths of several millimeters. 



Photoacoustic imaging has the potential to become a 
routinely used medical imaging modality, combining the 
superior contrast of optical techniques, with the penetra- 
tion depth of ultrasound.f^Its inherent ability to distin- 
guish regions of contrasting optical absorption make it an 
ideal candidate for imaging vascular structure, with pos- 
sible applications in diagnosis of stroke^ and early stage 
cancers 

Photoacoustic Tomography (PAT) is potentially ca- 
pable of producing real time 3D, high resolution im- 
ages to depths of several centimeters.^ The requirements 
of the ultrasonic detection system to achieve this are 
formidable, and as yet, no system has emerged which 
satisfy all criteria simultaneously. 

Though much emphasis has been placed on detector 
sensitivity, there are other, equally pressing requirements 
of a high resolution system. These include the need for 
a large detection surface, with high spatial resolution. 
Of particular importance for a system that can be used 
in clinical context is the ability to rapidly capture data, 
ideally over the whole detector surface simultaneously. 
A detector which is transparent to the excitation light 
is also favourable, as this allows a large amount of opti- 
cal energy to be dumped uniformly on the region being 
imaged. 

Piezoelectric detectors struggle with many of these 
requirements,^ and so a wide variety of optical detectors 
have been developed.^ Planar Fabry- Perot based sys- 
tems show good sensitivity and bandwidth response, but 
are typically slow to acquire data because of the need to 
tune the probing laser at each point on the detector sur- 
face to achieve peak sensitivity.^^ Simultaneous 2D data 
collection has been shown with these systems, but large 
detection areas are challenging to produce because of the 
difficulty in creating polymer coatings of uniform thick- 
ness. 

Microring resonators have been made with impressive 
sensitivity, element size and can be made transparent. 



However coupling and addressing a large array of micror- 
ings will be difficult, so it is yet to be seen if simultaneous 
full surface detection is achievable.^^ 

Pressure dependent optical reflectance detectors have 
been demonstrated with the ability to capture pressure 
data over a whole surface simultaneously,^^ without 
the need for complicated nanofabrication techniques of 
some other methods. Moreover, use of fast-framing, or 
gated CCDs simplifies data collection, and allows for high 
spatial resolution over a large detection surface. How- 
ever, the detection sensitivity of this type of system has 
so far been only modest, and may be insufficient for high 
resolution biological tissue imaging. 

We have developed a detector based on a Michelson in- 
terferometer (MI) with quadrature phase detection. This 
detector has comparable sensitivity to other optical de- 
tectors in the literature, but has the potential to per- 
form high resolution measurements over a full 2D sur- 
face simultaneously, without the need for any position 
dependent sensitivity tuning. The MI acts as an ultra- 
sound sensor simply by acoustically coupling an ultra- 
sound source to a mirror in one of the arms. The acoustic 
wave of pressure p, has an associated particle displace- 
ment ^, which shifts the position of the mirror as the 
wave passes through it. This change in position adjusts 
the relative phase of the laser beams in the two arms, re- 
sulting in a change in fringe brightness at the output of 
the interferometer. For a small amplitude wave travelling 
in the x direction at time pressure and displacement 
are related by: 



p{x,t) = -E 



dx ' 



(1) 



where E is the appropriate modulus of elasticity for the 
medium. 

The intensity, / of the recombined beam in a standard 
MI varies sinusoidally with mirror position: 



i = ^{i+cosm, 



(2) 
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where (j) = 47rn^/A. Here, Iq is the input intensity, (j) is 
the phase, A is the wavelength of the probe beam, and n 
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is the refractive index of the arm where the mirror posi- 
tion is changing. An MI therefore suffers from the same 
problem as Fabry- Perot type detectors, in that the laser 
(or mirror position) must be tuned to a sensitive region 
at each point in order to get good optical modulation for 
a given mirror displacement. The need for tuning was re- 
moved from our system by employing quadrature phase 
detection. In quadrature phase detection, two orthogo- 
nal linear polarizations are used to simultaneously obtain 
two separate interference patterns at the output of the 
interferometer, which have a relative phase difference of 
7r/2. This phase difference insures that the interference 
pattern of at least one of the polarization components is 
always in a sensitive region. 

The phase shift between the two polarization compo- 
nents is created by first linearly polarizing the light at 
45° from the vertical or horizontal axis of the polariz- 
ing beam splitter. A liquid crystal variable waveplate 
is placed in one of the arms, which retards the phase 
of one polarization component relative to the other by 
nominally 7r/4 in both the forward and reverse trips. A 
variable retarder is used instead of a fixed A/8 plate to 
compensate for small amounts of birefringence present in 
other optical components. 

The 'phase sensitivity', |^ (which is the optical in- 
tensity modulation per radian of phase) of a normal MI 
varies between rad~^ and Iq/2 rad~^. For an MI with 
quadrature detection, the phase sensitivity for each po- 
larization is simply added together, so is always between 
/o/2 rad~^ and Iq/\/2 rad~^. This ensures that the to- 
tal sensitivity of the system is always at least as high as 
the maximum of a standard MI, irrespective of absolute 
mirror position. 

A diagram of the detector setup can be seen in Fig. 
[l] The mirror position is recovered from the detected 
intensities of the two polarizations, /i^2 by first scaling 
them between —1 and 1, then treating them as points on 
the unit circle: (j) = atan2{Ii^ I2). 

Our setup uses an expanded 5 mW CW 633 nm HeNe 
laser as the probe, and the detectors are amplified pho- 
todiodes with a 20 MHz bandwidth. The photodiodes 
are apertured, which sets the spatial resolution of the 
system. The signal of each photodiode is recorded using 
a digital oscilloscope. The sensing mirror is a 150 /im 
thick glass substrate with a gold reflective coating. To 
maintain mechanical stability, the mirror is mounted on 
an optical window by bonding it around the perimeter 
with epoxy resin. This method of bonding also insures 
the presence of an air gap between the reflective surface 
and the window, eliminating the possibility of acoustic 
waves propagating into the window. The glass surface of 
the mirror is acoustically coupled to the medium being 
imaged, with the window becoming part of the interfer- 
ometer arm. The resulting mirror can be seen in Fig. 

m 

It is necessary to use such a thin mirror for two reasons. 
Firstly, a thick mirror suffers from the acoustic wave re- 
flecting back and forth off the boundary of the substrate. 
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FIG. 1. A schematic of the Michelson interferometer detec- 
tor. Ultrasound passing through the sensing mirror changes 
the path length of the light, which alters the phase of the 
interference pattern. Quadrature phase detection allows high 
sensitivity to be achieved, whatever the mirror positions. 
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FIG. 2. (a) Schematic of the sensing mirror showing probe 
laser, (b) The back side of the sensing mirror used in this 
experiment. 



This reflected wave interferes with the incoming wave, 
making the detected displacement useless. Multiple re- 
flections are strongly suppressed in a mirror that is much 
thinner than the acoustic wavelength, so the displace- 
ment on the mirror surface accurately represents only 
the incoming wave. Secondly, a mirror of thickness sig- 
nificantly greater than the acoustic wavelength is able to 
support surface waves. These surface, or Rayleigh waves 
can exist whenever there is an impedance mismatch be- 
tween two media, and are generated on the surface of 
the mirror when the photoacoustic pulse first reaches the 
boundary. The Rayleigh waves then propagate outwards 
across the surface and interfere with incoming photoa- 
coustic waves. 

To see the fundamental limitations of the detection 
system, it is useful to look at the theoretically achiev- 
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able sensitivity. Like all Fabry-Perot and piezoelectric 
type detectors, the MI detector is sensitive to particle 
displacements, rather than directly to pressure. However 
because the MI detects the absolute position of a single 
plane, rather than the relative position of two planes, its 
sensitivity to pressure is easier to describe analytically. 
For a sinusoidal acoustic wave of pressure amplitude po^ 
travelling in the positive x direction, the corresponding 
particle displacement is given by: 



Po 

27TUZ 



cos(kx — 27rut), 



(3) 



where k is wavenumber, u is frequency, and z is the spe- 
cific acoustic impedance of the propagation medium. For 
optical based pressure sensors the sensitivity, S can be 
simply given as a proportion of optical intensity modu- 
lation per unit acoustic pressure: S = I/Iq^. Using the 

dcj) dp' 

is always at least 0.5/o, and 



chain rule, this may be expanded to ^ 
the quadrature MI, 



dl 



d^ 



is simply calculated from the expression for (p. g can 
be calculated from Eqn. [3| however it must be modified 
to describe the setup employed in our system. Firstly, 
since the pressure wave must propagate from the orig- 
inal medium into the glass mirror substrate, the pres- 
sure must be multiplied by the transmission coefficient: 
Tp = 2z2/{zi + Z2), where Zi^2 are the specific acoustic 
impedances of the first and second media respectively.'^ 
Also, the mirror is essentially on a free boundary (since 
Zair Zgiass)i SO the particle displacement will be twice 
as great as in the bulk. Combining these terms gives the 
expression for the frequency dependent sensitivity: 



An 



vX{zi + Z2) ' 



(4) 



Typical values of acoustic impedances for water and glass 
are 1.5 x 10^ Pa s m~^ and 13.1 x 10^ Pa s m~^ respec- 
tively. Taking n = 1 as the refractive index of air, and 
letting A = 633 nm, the sensitivity of our detector is 
S{v) = Hz Pa~^. For a photoacoustic wave of 

frequency 1 MHz, the sensitivity of the detector should 
be 4.3% optical modulation per 100 kPa of peak acous- 
tic pressure. This value can be compared directly with 
the sensitivity for an optical refiectance based detector 
of 0.19 to 0.81% reported elsewhere in the liter at ure.'I^ 

Fig. [3] shows the actual detected displacement caused 
by a PA wave from a single source, positioned 3.3 mm 
directly behind the detector. To recover the acoustic 
pressure from the displacement, the temporal derivative 
must be taken according to Eqn. [l] where c = dx/dt 
has been used to change the variable of differentiation. 
A post processing low pass filter is applied to the de- 
tected displacement before the derivative is taken. This 
ensures the calculated pressure does not contain unphys- 
ical spikes which are artefacts of taking the derivative of 
a noisy signal. 

The noise equivalent displacement in this trace is ap- 
proximately 5 A, where the signal has been averaged for 
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FIG. 3. (a) The detected mirror displacement due to a pho- 
toacoustic wave produced from a single source and (b), the 
corresponding pressure wave. A comparison between simu- 
lated and experimental pressure data shows the calculated 
sensitivity for the system is accurate. 



64 pulses. The majority of the noise in this case is due to 
the amplified photodiodes, which may be improved using 
different equipment. 

The comparison between a simulation and an actual 
detected photoacoustic wave are shown in Fig. [3) and the 
agreement between the two support the values of sensi- 
tivity calculated previously. In both the simulation and 
experiment, a 500 /im diameter tube was filled with di- 
luted India ink with absorption coefficient 35 cm~^ at 
1064 nm, which is similar to the optical absorption coef- 
ficient of blood at wavelengths commonly used in PAT. 
The tube was illuminated with a 10 ns pulse of 1064 nm 
light with ffuence of 25 mJ cm~^. The resulting PA wave 
was detected directly in front of the cylinder. The simu- 
lation assumed an acoustically homogeneous water prop- 
agation medium. For the experimental measurement, the 
tube was submerged in a water filled glass cell. One wall 
of the cell was made of acoustically transparent polyethy- 
lene film, which was coupled to the sensing mirror of the 
interferometer with commercial ultrasound coupling gel. 
The detected pressure represents the pressure in the bulk 
of the glass substrate. This value was divided by the pres- 
sure transmission coefficient to give the acoustic pressure 
in the water, allowing a direct comparison with simula- 
tion. 

To create photoacoustic images, the photodiodes could 
be moved laterally to build up ID or 2D scans. However 
the lower laser intensity at the edges of the expanded 
probe beam would mean that data collected in these re- 
gions had a lower signal to noise ratio (SNR). Instead, 
the sample itself was scanned laterally. This allowed dif- 
ferent regions of the generated ultrasonic wavefield to be 
sampled in an equivalent manner to scanning the photo- 
diodes, while maintaining a high SNR. A ID scan of the 
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FIG. 4. (a) Detected acoustic pressure from a single photoa- 
coustic source and (b) , the corresponding reconstructed image 
of the source. 



wavefield produced by the cyHndrical optical absorber is 
shown in Fig. [4j The source was diluted India ink in a 
500 /am internal diameter silicone tube, and was illumi- 
nated by the same pulsed laser as previously described. 
The photodiode aperture was set to a diameter of 200 /im, 
which was the same as the lateral step size. 

The reconstruction algorithm used here was described 
by Kostli et al^ and implemented in MATLAB® us- 
ing the k-Wave photoacoustic package.^ It is based on 
Fourier transforms, and it is theoretically exact if pres- 
sure is detected over an infinite plane for infinite time. In 
the reconstruction shown in Fig. [4j the source has been 
positioned correctly, but has suffered some blurring and 
distortion which is consistent with other implementations 
of this inversion technique. The blurring is unsurprising 
given that the diameter of the source was only 500 jam 
and the spatial separation of each data point was 200 /im. 

The resolution of our system is currently limited by the 
size of the aperture in front of the photodiodes, which in 
turn is limited by the need to get sufficient laser power 
to the photodiodes. This is easily improved by increas- 
ing the laser power, and could be achieved using inex- 
pensive diode lasers. Using a shorter wavelength probe 
laser would also be a simple way to boost sensitivity ac- 
cording to Eqn. [4j Manually scanning the photodetec- 
tors (or the source) to build up an image is too slow 
for real-time imaging applications, so any useful system 
must ultimately be capable of performing simultaneous 
detection over the whole surface. This could be simply 
achieved in our system by replacing the photodiodes with 
fast gated intensified charge-coupled device (ICCD) cam- 
eras. The potential to use CCDs is a significant advan- 



tage over other proposed optical photoacoustic detectors, 
which have no easy route to simultaneous measurement 
of all elements over a large surface. 

The current configuration of our sensing mirror (con- 
structed from a thin glass substrate) may cause limi- 
tations to the detectable bandwidth needed in higher 
resolution systems, due to the possible reappearance of 
Rayleigh waves at higher frequencies. However this could 
be addressed by using polymer substrates impedance 
matched to water. Also, where deeper imaging is re- 
quired, the bandwidth requirements of the detector are 
much more forgiving, since very high frequency acous- 
tic waves are strongly attenuated in tissue. As such, the 
current sensing mirror is suitable for imaging depths be- 
yond a few millimeters with no reduction in attainable 
resolution. 

In summary, we have demonstrated the use of a Michel- 
son interferometer as a photoacoustic detector. The use 
of quadrature phase detection removes the need to tune 
the sensitivity of the interferometer, allowing for the 
possibility of simultaneous full surface detection. We 
used the detector to produce proof-of-concept photoa- 
coustic images with sub-millimeter resolution and sug- 
gested ways this could be improved. Future work will in- 
volve demonstrating simultaneous full surface measure- 
ment, with the aim of producing real time 3D photoa- 
coustic visualizations. 
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